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Abstract 

A search is performed for narrow resonances decaying into WW, WZ, or ZZ boson pairs 
using 20.3 fb“* of proton-proton collision data at a centre-of-mass energy of = 8 TeV 
recorded with the ATLAS detector at the Large Hadron Collider. Diboson resonances with 
masses in the range from 1.3 to 3.0 TeV are sought after using the invariant mass distribution 
of dijets where both jets are tagged as a boson jet, compatible with a highly boosted W or Z 
boson decaying to quarks, using jet mass and substructure properties. The largest deviation 
from a smoothly falling background in the observed dijet invariant mass distribution occurs 
around 2 TeV in the WZ channel, with a global significance of 2.5 standard deviations. 
Exclusion limits at the 95% confidence level are set on the production cross section times 
branching ratio for the WZ final state of a new heavy gauge boson, W', and for the WW 
and ZZ final sfafes of Kaluza-Klein excifafions of fhe gravifon in a bulk Randall-Sundrum 
model, as a funcfion of tbe resonance mass. W' bosons with couplings predicted by the 
extended gauge model in the mass range from 1.3 to 1.5 TeV are excluded at 95% confidence 
level. 
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1 Introduction 


The substantial dataset of Large Hadron Collider (LHC) proton-proton (pp) collisions at = 8 TeV 
collected by the ATLAS experiment provides a distinct opportunity to search for new heavy resonances 
at the TeV mass scale. This paper presents a search for narrow diboson resonances {WW, fTZ and ZZ) 
decaying to fully hadronic final states. The fully hadronic mode has a higher branching fraction than 
leptonic and semileptonic decay modes, and is therefore used to extend the reach of the search to the 
highest possible resonance masses. 

W and Z bosons resulting from the decay of very massive resonances are highly boosted, so that each 
boson’s hadronic decay products are reconstructed as a single jet. The signature of the heavy resonance 
decay is thus a resonance structure in the dijet invariant mass spectrum. The dominant background for this 
search is due to dijet events from QCD processes, which produce a smoothly falling spectrum without res¬ 
onance structures. To cope with this large background, jets are selected using a boson tagging procedure 
based on a reclustering-mass-drop filter (BDRS-A, similar to the method introduced in ref. [1]), jet mass, 
and further substructure properties. The tagging procedure strongly suppresses the dijet background, al¬ 
though these QCD processes still overwhelm the expected backgrounds from single boson production 
with one or more jets. Standard Model (SM) diboson production, single-top and top-pair production. As 
all of these background sources produce dijet invariant mass distributions without resonance peaks, the 
expected background in the search is modelled by a fit to a smoothly falling distribution. 

Diboson resonances are predicted in several extensions to the SM, such as technicolour [2-4] , warped 
extra dimensions [5-7], and Grand Unified Theories [8-11]. To assess the sensitivity of the search, to 
optimise the event selection, and for comparison with data, two specific benchmark models are used: an 
extended gauge model (EGM) W —> WZ where the spin-1 W gauge boson has a modified coupling fo 
fhe SM W and Z bosons [12-14], and a spin-2 graviton, Grs WW or ZZ, a Kaluza-Klein mode [5, 15] 
of fhe bulk Randall-Sundrum (RS) graviton [16-18]. 

The CMS collaborafion has performed a search for diboson resonances wifh fhe fully hadronic final 
sfafe [19] of comparable sensifivify to fhe one presenfed in fhis article. In fhis search, fhe EGM W' —> WZ 
wifh masses below 1.7 TeV and Grs of fhe original RS model decaying fo WW wifh masses below 1.2 
TeV are excluded af 95% confidence level (CE). The CMS collaborafion has also published upper limits 
on the production of generic diboson resonances using semileptonic final sfafes [20]. Using fhe final 
sfafe, fhe ATEAS collaborafion has excluded af 95 % CE a bulk Grs —> ZZ wifh mass below 740 GeV [21]. 
Eor narrow resonances decaying exclusively fo WZ or WW, fhe sensifivify of fhe ATEAS search in fhe 
£vqq channel [22] is comparable fo fhaf of fhe search presenfed here in fhe mass range from 1.3 fo 2.5 
TeV. Thai search has also excluded a bulk Grs —> WW wifh mass below 760 GeV. 

2 ATLAS detector and data sample 

The ATEAS defecfor [23] surrounds nearly fhe enfire solid angle around fhe ATEAS collision poinf. 
If has an approximately cylindrical geomefry and consisfs of an inner fracking defecfor surrounded by 
elecfromagnelic and hadronic calorimelers and a muon speclromeler. The fracking defecfor is placed 
wifhin a 2 T axial magnefic field provided by a superconducfing solenoid and measures charged-parficle 
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trajectories with pixel and silicon microstrip detectors that cover the pseudorapidity ^ range \t]\ < 2.5, and 
with a straw-tube transition radiation tracker covering \r]\ < 2 . 0 . 

A high-granularity electromagnetic and hadronic calorimeter system measures the jets in this analysis. 
The electromagnetic calorimeter is a liquid-argon (LAr) sampling calorimeter with lead absorbers, span¬ 
ning I 77 I <3.2 with barrel and end-cap sections. The three-layer central hadronic calorimeter comprises 
scintillator tiles with steel absorbers and extends to I 77 I = 1.7. The hadronic end-cap calorimeters measure 
particles in the region 1.5 < \t]\ < 3.2 using liquid argon with copper absorber. The forward calori¬ 
meters cover 3.1 < \q\ <4.9, using LAr/copper modules for electromagnetic energy measurements and 
LAr/tungsten modules to measure hadronic energy. 

Events are recorded in ATLAS if they satisfy a three-level trigger requirement. The level-1 trigger detects 
jet and particle signatures in the calorimeter and muon systems with a fixed latency of 2.5 /rs, and is 
designed to reduce the event rate to less than 75 kHz. Jets are identified af level-1 wifh a sliding-window 
algorifhm, searching for local maxima in square regions wifh size Arjx = 0.8 x 0.8. The subsequenf 
high-level frigger consisfs of fwo sfages of soflware-based trigger filters which reduce the event rate to a 
few hundred Hz. Events used in this search satisfy a single-jet trigger requirement, based on at least one 
jet reconstructed at each trigger level. At the first filtering stage of the high-level trigger, jet candidates are 
reconstructed from calorimeter cells using a cone algorithm with small radius, /?=0.4. The final filler in 
Ihe high-level frigger requires a jel lo safisfy a higher Iransverse momenlum (pj) Ihreshold, reconslrucled 
wifh Ihe anli-k/ algorifhm [24] and a large radius paramefer (7?=1.0). 

This search is performed using Ihe dafasef collecfed in 2012 from 8 TeV LHC pp collisions using a 
single-jel frigger wifh a nominal pj Ihreshold of 360 GeV. The inlegraled luminosily of Ihis dafasef afler 
requiring good beam and defector conditions is 20.3 fb“^, wifh a relative uncerfainly of +2.8%. The 
uncerfainfy is derived following Ihe mefhodology defailed in ref. [25]. 


3 Simulated data samples 

The leading-order Monle Carlo (MC) generator Pythia 8.170 [26] is used to model W WZ evenls in 
order lo defermine and opfimise fhe sensilivify of fhis search. Pythia 8 uses fhe px-ordered showering 
inlroduced in Pythia 6.3 [27, 28], and inlerleaves mulfiple parlon inleraclions wifh bolh inilial- and final- 
slale radiation. The samples generaled for Ibis analysis use MSTW2008 [29] parlon dislribulion funcfions 
(PDEs), wifh parlon shower parameters luned to ATLAS underlying-evenl dala [30]. Hadronisalion is 
based on fhe Lund siring fragmenfalion framework [31]. An addilional sef of W samples generaled wifh 
Pythia for fhe hard scattering interaction and Herwig++ [32] for parlon showering and hadronisalion is 
used fo assess syslemalic uncerlainfies on fhe signal efficiency due to uncerlainfies on fhe parlon shower 
and hadronisalion model. These samples use angular-ordered showering and clusfer hadronisalion. 

The W boson samples are generaled for differenl resonance masses, covering fhe range 1.3 < mw < 3.0 
TeV in 100 GeV inlervals. The W is required lo decay lo a IT and a Z boson, which are bolh forced lo 
decay hadronically. The cross secfion times branching ralio as well as fhe resonance widlh for fhe samples 
lisfed in fable 1 are calculafed by Pythia 8 assuming EGM couplings [12] for fhe W. In parlicular, fhe W 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the i/-axis points 
upward. Cylindrical coordinates (r, 0) are used in the transverse plane, (p being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 0 as ?/ = - In tan(0/2). 
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coupling to VKZ is equal to that of the W coupling scaled by c x where c is a coupling scaling 

factor of order one which is set to unity for the samples generated here. The partial width of W WZ 
decays thus scales linearly with mvrs leading to a narrow width over the entire accessible mass range. 
Because of the anti-quark parton distribution functions involved in the production, a significant part of 
the W cross section for large W masses is due to off-shell inferacfions which produce a low-mass fail in 
fhe W mass specfrum. The relafive size of fhe low-mass fail increases wifh fhe W mass: fhe fraction of 
evenfs wifh a diboson mass below 20% of fhe pole mass of fhe W increases from 10% for mw = 1.3 TeV 
fo 22% for myfi = 2.0 TeV and fo 65% for mw' = 3.0 TeV. 

An exfended RS model wifh a warped exfra dimension is used for fhe excited graviton benchmarks. 
In Ibis model fhe SM fields are allowed fo propagate in fhe warped exfra dimension [16], avoiding fhe 
consfrainfs on fhe original RS model from limifs on flavour-changing neufral currenfs and elecfroweak 
precision measuremenfs. The model is characterised by a dimensionless coupling consfanf k/Mpi ~ 1, 
where k is fhe curvafure of fhe warped exfra dimension and Mpi is fhe reduced Planck mass. The RS 
excifed gravifon samples are generated wifh CalcHEP 3.4 [33] selling kjMpi = R covering fhe resonance 
mass range 1.3 < nicRs < 3.0 TeV in 100 GeV intervals. The graviton resonance is decayed fo WW or 
ZZ, and fhe resulfing IT or Z bosons are forced to decay hadronically. The cross seclion limes branching 
rafio as well as fhe resonance widfh calculaled by CalcHEP for fhe RS model are lisfed in fable 1. Evenfs 
are generaled using CTEQ6E1 [34] PDEs, and use Pythia 8 for fhe parlon shower and hadronisafion. 

To characferise fhe expected dijel invarianf mass specfrum in fhe mass range 1.3-3.0 TeV, simulafed QCD 
dijel evenfs, diboson evenfs, and single IT or Z bosons produced wifh jefs are used. Confribulions from 
SM diboson evenfs are expecled to accounl for approximately 6% of fhe selecled sample, and single 
boson produclion is expecled to conlribule less lhan 2%. Confribulions from ti producfion, sludied us¬ 
ing MC@NEO [35] and HERWIG [36] showering, were found fo be negligible and are nol considered 
furlher. 

QCD dijel evenfs are produced wifh Pythia 8 and fhe CTIO [37] PDEs and fhe W/Z -i- jefs samples are 
produced wifh Pythia 8 and CTEQ6E1 PDEs. Diboson evenfs are produced al fhe generator level wifh 
POWHEG [38], using Pythia for fhe sofl parlon shower. The samples of single IT or Z bosons produced 
wifh jefs are furlher used to delermine a scale faclor for fhe efficiency of fhe boson fagging seleclion, by 
comparing fhe boson fagging efficiencies belween simulafion and collision dafa in a IT/Z-i-jels-enriched 
sample. 

The final-sfale particles produced by fhe generators are propagaled fhrough a defailed defector simula- 
lion [39] based on GEANT4 [40]. The average number of pp interactions per bunch crossing was ap- 
proximafely 20 while fhe collision dafa were collected. The expected confribulion from fhese additional 
minimum-bias pp inferacfions is accounled for by overlaying additional minimum-bias evenfs generated 
wifh Pythia 8 , malching fhe dislribufion of fhe number of inferacfions per bunch crossing observed in 
collision dafa. Simulafed evenfs are fhen reconstrucled wifh fhe same algorifhms run on collision dafa. 


4 Boson jet identification 


In Ihis search IT and Z bosons from fhe decay of fhe massive resonance are produced wifh a large frans- 
verse momenlum relafive to Iheir mass and each boson is reconslrucled as a single large-radius jel. Boson- 
jef candidates are fhen idenfified by applying lagging requiremenls based on fhe reconslrucled jel proper- 
lies, as described below. 
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Table 1: The resonance width (T) and the product of cross sections and branching ratios (BR) to four-quark final 
states used in modelling W' —» WZ, Grs ^ WW, and Grs —> ZZ, for several values of resonance pole masses (m). 
The fraction of events in which the invariant mass of the W or Grs decay products lies within 10% of the nominal 
resonance mass (/io%) is also displayed. 


m 

[TeV] 

Yw 

[GeV] 

^"Grs 

[GeV] 

IT' ^ 
cr xBR 
[fb] 

ITZ 

/io% 

Grs ^ 
cr xBR 
[fb] 

WW 

fio% 

Grs ^ 
cr xBR 
[fb] 

ZZ 

fio% 

1.3 

47 

76 

19.1 

0.83 

0.73 

0.85 

0.37 

0.84 

1.6 

58 

96 

6.04 

0.79 

0.14 

0.83 

0.071 

0.84 

2.0 

72 

123 

1.50 

0.72 

0.022 

0.83 

0.010 

0.82 

2.5 

91 

155 

0.31 

0.54 

0.0025 

0.78 

0.0011 

0.78 

3.0 

109 

187 

0.088 

0.31 

0.00034 

0.72 

0.00017 

0.71 


4.1 Jet reconstruction 

Jets are formed by combining topological clusters [41] reconstructed in the calorimeter system, which are 
calibrated in energy with the local calibration scheme [42] and are considered massless. These topological 
clusters are combined into jets using the Cambridge-Aachen (C/A) algorithm [43, 44] implemented in 
FastJet [45] with a radius parameter R = 1.2 with a radius parameter R = 1.2. The C/A algorithm 
iteratively replaces the nearest pair of elements (topological clusters or their combination) with their 
combination^ until all remaining pairs are separated by more than R, defining the distance A/? between 
elements as (A/?)^ = (A^)^ -i- (A0)^ where y is the rapidity. The jets are the elements remaining after this 
final sfage of iferafion, and fhe lasf pair of elemenfs fo be combined info a given jef are referred fo here as 
fhe subjefs of fhaf jef. Charged-parficle fracks reconstrucfed in fhe fracking defector are mafched fo fhese 
calorimefer jefs if fhey fall wifhin fhe passive cafchmenf area of fhe jef [46], determined by representing 
each frack by a collinear “ghosf” consfifuenl wifh negligible energy during jef reconsfrucfion. Only well- 
reconsfrucfed fracks wifh pj > 500 MeV and consisfenf wifh parficles originafing af fhe primary^ collision 
verfex are considered. 

The jefs are fhen groomed fo idenfify fhe pair of subjefs associafed wifh fhe W qq' or Z ^ qq decay, 
and fo reduce fhe effecf of pileup and ofher noise sources on fhe resolufion. The grooming algorifhm is a 
varianf of fhe mass-drop filtering fechnique [ 1 ], which firsl examines the sequence of pairwise combina¬ 
tions used to reconstruct the jet in reverse order. At each step the lower-mass subjet is discarded, and the 
higher-mass subjet is considered as the jet, continuing until a pair is found which satisfies mass-drop and 
subjet momentum balance criteria parameterised by pf and respectively. Iteration stops when a pair 
of subjets is found for which each subjet mass mq-) satisfies p = my^lmo < pf, and for which 

ARn^ b) 

^Jy = minipTj pTj ) -^ > Vw , 

where pjj^ and pjj^ are the transverse momenta of subjets ji and 72 respectively, AR(j^Jp is the distance 
between subjets ji and 72 , and mo is the mass of the parent jet. 

^ Elements are combined by summing their four-momenta. 

^ The primary collision vertex is the reconstructed vertex with the greatest sum of associated track 
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The subjet momentum-balance threshold y{ and the mass-drop parameter /if used in this analysis, given 
in table 2, are chosen to stop the iteration when the two subjets corresponding to the W or Z boson decay 
have been identified in simulated signal events. The best signal to background ratio for a given signal 
efficiency was obfained wifh /if sef fo 1, hence no mass-drop requiremenf is applied in fhis analysis. 

The selected pair of subjefs is fhen filfered: fhe original fopological cluster consfifuenfs of fhaf pair of 
subjefs are faken fogefher and clusfered using fhe C/A algorifhm wifh a small radius parameter (R^ = 0.3), 
and all buf fhe fhree = 3) highesf-pj jets resulting from this reclustering of the subjets’ constituents 
are discarded. If there are fewer than three jets after the reclustering, all constituents are kept. Those 
constituents that remain form the resulting filtered jet, which is further calibrated using energy- and ly- 
dependent correction factors derived from simulation by applying a procedure similar to the one used 
in ref. [47]. The calibrated four-momentum is used as the IT or Z boson candidate’s four-momentum in 
subsequent cuts and in reconstructing the heavy resonance candidate’s mass in each selected event. 

Table 2; Parameters for the mass-drop filtering algorithm used to groom C/A jets. The choice of /if parameter 
corresponds to no mass-drop requirement being imposed in the grooming procedure. 


Filtering parameter 

Value 


0.2 

Rt 

1 

Rr 

0.3 

nr 

3 


4.2 Boson jet tagging 

The grooming algorithm rejects jets that do not satisfy the momentum balance and mass-drop criteria 
at any stage of iteration, and thus provides a small degree of discriminating power between jets from 
hadronically decaying bosons and those from QCD dijet production. To improve the discrimination 
in this analysis, the remaining jets are also tagged with three additional boson tagging requirements. 
First, a more stringent subjet momentum-balance criterion {^/y > 0.45) is applied to the pair of subjets 
identified by the filtering algorithm at the stopping point before the reclustering stage, since jets in QCD 
dijet events that survive grooming tend to have unbalanced subjet momenta characteristic of soft gluon 
radiation. Figure 1(a) shows the subjet momentum-balance distribution for jets in signal and QCD dijet 
background simulated events. Unlike jets from massive boson decays in which the hadron multiplicity 
is essentially independent of the jet pj, energetic gluon jets are typically composed of more hadrons, 
so the number of charged-particle tracks associated with the original, ungroomed jet is required to be 
small (iifrk < 30). Figure 1(b) shows the number of tracks matched to jets for selected jets in signal and 
background simulated events. The efficiency of this selection requirement must be corrected by a scale 
factor derived from data, as explained in section 7. Finally, a selection window is applied to the invariant 
mass of the filtered jet, mj, since this quantity is expected to be small for jets in QCD dijet events and to 
reflect the boson mass for jets from hadronic boson decays. The expected mass distribution of jets in Grs 
simulated events and the dijet background simulation is illustrated in figure 1(c). Narrow mass windows 
with a width of 26 GeV are chosen to optimise sensitivity to signal events and are centred at either 82.4 
GeV or 92.8 GeV, where the mass distributions of the W and Z jets, respectively, peak in simulation. Each 
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jet’s mass must fall within either the IV or Z mass window, consistent with the IVZ, WW or ZZ final state 
being studied. 



Jet '{y 



(a) 


(b) 



(c) 

Figure 1; The distribution of the boson-tagging variables (a) subjet momentum balance sjy, (b) number of tracks 
Htrk matched to the jet, and (c) mass mj of the groomed jet, in simulated signal and background events. The signal 
and background distributions are normalised to unit area, and the last bin of each histogram includes the fraction 
of events falling outside of the displayed range. Requirements are placed on the events used to ensure that the 
kinematics of the signal and background events are comparable. 


5 Event selection 

High-mass resonances decaying to a pair of boosted vector bosons with subsequent hadronic decay are 
recognised as two large-radius massive jets with large momentum, typically balanced in pj. Events 
in this search must therefore first satisfy the high-pj large-radius jet trigger, which is found to select 
over 99% of C/A R-1.2 jets within \t]\ < 2.0 and with ungroomed pj greater than 540 GeV. Events 
are removed if they contain a prompt electron candidate with Ej > 20 GeV in the regions \t]\ < 1.37 or 
1.52 < It/I < 2.47, or a prompt muon candidate with px > 20 GeV in the region \r]\ < 2.5. This requirement 
ensures that this analysis has no events in common with other diboson search analyses [21, 22]. Events 
with reconstructed missing transverse momentum exceeding 350 GeV are also removed, as these are used 
in searches sensitive to diboson resonances with a Z boson decaying to neutrinos [48] . 


































5.1 Event topology 


For events satisfying the requirements above, two C/A R=\.l jets with pj exceeding 20 GeV must be 
found and must pass the mass-drop filtering procedure. The two jets with the highest transverse mo¬ 
mentum must have |? 7 | < 2.0 to ensure sufficient overlap with the inner tracking detector, since associated 
charged-particle tracks are used in the boson tagging requirements and for estimating systematic uncer¬ 
tainties. In addition, a requirement on the rapidity difference between the two leading jets, \y\ - y 2 \ < 1-2, 
is imposed to improve the sensitivity. This rapidity difference is smaller for ^-channel processes such as 
the W and Grs signal models than for the t-channel processes dominating the QCD dijet background. 

The combined efficiency of these three cuts in W signal events is between 72% and 81% depending on 
the resonance mass, for events from each signal sample in which the true diboson mass lies within 10% 
of the nominal value. For Grs signal events, the combined efficiency is between 82% and 87%. The 
difference in the expected efficiencies between the W and Grs signals is related to the different event 
topologies expected for spin-1 and spin-2 resonances affecting acceptance. 

A selection on the pj asymmetry of the two leading jets, (pxi - PT 2 )liPTi + Pti) < 0.15, is used to 
reject events where one of the jets is poorly measured or does not come from the primary pp collision. 
The signal selection efficiency of this cut exceeds 97% in W signal samples and 90% in the samples 
of Grs events. The difference in the expected efficiencies between the W and Grs signals is related 
to their different production mechanisms. Figure 2(a) shows the selection efficiency of the event topo¬ 
logy requirements for signal events with resonance mass within 10% of the nominal signal mass for the 
W WZ, bulk Grs — > WW and bulk Grs ^ ZZ benchmark models, with statistical and systematic 
uncertainties indicated by the width of the bands in the figure. 


ATLAS Simulation 
= 8 TeV 


□ egm W'^WZ 
^bulkG„3^WW 
I I bulk Gpjg —> ZZ 



1.4 1.6 


2.2 2.4 2.6 2.8 3 

Resonance Mass [TeV] 



(a) (b) 

Figure 2: Event selection efficiencies as a function of the resonance masses for EGM W WZ and bulk Grs ^ 
WW and ZZ for simulated events with resonance mass within 10% of the nominal signal mass. In (a), the event 
topology requirements are applied to EGM W' WZ, Grs ^ WW and Grs —» ZZ samples, while in (b), the WZ, 
WW and ZZ boson tagging selections are also applied in the EGM W' —> WZ, Grs ^ WW and Grs —> ZZ samples 
respectively and the efficiencies shown are corrected by the simulation-to-data scale factor. The width of the bands 
in each figure indicates both the statistical and systematic uncertainties. 
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5.2 Boson tagging requirements 

The two jets with the highest transverse momentum each must satisfy the three boson tagging require¬ 
ments discussed in section A'. > 0.45, ntrk < 30, and \mj - mv\ <13 GeV, where my is the peak 

value of the reconstructed IT or Z boson mass distribution. For the W ITZ search, this final cut sets 
my equal to the peak reconstructed W boson mass when applied to the lower mass jet, and to the peak 
reconstructed Z boson mass when applied to the higher mass jet. 

The expected efficiency of these boson tagging cuts applied to signal events is evaluated using the MC 
signal samples described in section 3. For signal events passing event topology requirements on the mass- 
drop filtering, r], the rapidity difference, and the py asymmetry, the average efficiency of the tagging cuts 
for each of the two leading-px filtered jets is approximately the same in the Grs ^ WW and Grs —> ZZ 
samples, and ranges from 44.0% in the wtcrs - 1-2 TeV sample to 33.9% for the /wcrs = 3.0 TeV sample. 
Figure 2(b) shows the selection efficiency of the event selection and tagging requirements for signal events 
with resonance mass within 10% of the nominal signal mass for the W —> WZ and bulk Grs —> WW 
and bulk Grs ^ ZZ benchmark models, with both statistical and systematic uncertainties included in 
the error band. The average background selection efficiency of the tagger for each of the two leading-px 
filtered jets in simulated QCD dijet events satisfying the same event selection requirements ranges from 
1.2% for events with dijet masses between 1.08 TeV and 1.32 TeV, to 0.6% for events with dijet masses 
between 2.7 TeV and 3.3 TeV. 

5.3 Dijet mass requirement 

The invariant mass calculated from the two leading jets must exceed 1.05 TeV. This requirement restricts 
the analysis of the dijet mass distribution to regions where the trigger is fully efficient for boson-tagged 
jets, so that the trigger efficiency does not affect its shape. 

6 Background model 

The search for high-mass diboson resonances is carried out by looking for resonance structures on a 
smoothly falling dijet invariant mass spectrum, empirically characterised by the function 


ax 


( 1 ) 


where x = nr/// and mjj is the dijet invariant mass, pi is a normalisation factor, p 2 and p^ are 
dimensionless shape parameters, and ^ is a dimensionless constant chosen after fitting to minimise the 
correlations between p 2 and p^. A maximum-likelihood fit, with parameters p\, P 2 and p^ free to float, is 
performed in the range 1.05 TeV < rujj <3.55 TeV, where the lower limit is dictated by the point where 
the trigger is fully efficient for tagged jets and the upper limit is set to be in a region where the data and 
the background estimated by the fit are well below one event per bin for the tagged distributions. The 
likelihood is defined in terms of events binned in 100-GeV-wide bins in rujj as 



( 2 ) 
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where n, is the number of events observed in the mjj bin and d,- is the background expectation for the 
same bin. 

The functional form in eq. 1 is tested for compatibility with distributions similar to the expected back¬ 
ground by applying it to simulated background events and to several sidebands in the data. Figure 3 
shows fits to the Herwig-i--i- and Pythia simulated dijet events that pass the full event selection and tag¬ 
ging requirements on both jets, where the predictions from these leading-order generators are corrected 
by reweighting the untagged leading-jet pj distributions to match the untagged distribution in data. Fig¬ 
ure 4 shows the results of fitting the dijet mass spectrum before tagging, and in otherwise tagged events 
where both the leading and subleading jet have masses falling below the boson-tagging mass windows, 
in the range 40 < mj < 60 GeV. For the data selected before boson tagging, the trigger efficiency as a 
function of the dijet mass is taken into account in the fit, because in the untagged jet sample the trigger 
is not fully efficient in the first dijet mass bin displayed. The fitted background functions in figures 3 
and 4 are integrated over the same bins used to display the data, and labelled “background model” in the 
figures. The size of the shaded band reflects the uncertainties of the fit parameters. In figure 4(a) the 
background model is shown, but the size of the uncertainty band is too small to be seen. The lower insets 
in the figures show the significance, defined as the signed z-value of the difference befween the distri¬ 
bution being modelled and the background model’s prediction [49]. The significance with respect to the 
maximum-likelihood expectation is displayed in red, and the significance when taking the uncertainties 
on the fit parameters into account is shown in blue. 

Table 3 summarises the results of these fits, as well as fits to data where one jet mass falls in the low-mass 
sideband (40 < mj < 60 GeV) and the other falls in a high-mass sideband from 110 < mj < 140 GeV, 
and where both jet masses fall in the high-mass sideband. The dijet mass distribution of the simulated 
background and of each of these background-dominated selections are well-described by the functional 
form in eq. (1). 


Table 3: Goodness-of-fit for maximum-likelihood fits of the background model to the dijet mass distribution in 
simulated events, and in selected mass sidebands from data events where at least one of the leading and subleading 
jet fails the jet mass selection. One-sided probabilities are displayed; for the three data sideband fits, these 
probabilities were calibrated using pseudo-experiments to avoid biases due to empty bins. 


Sample 

vVnDOF 

Probability 

Pythia dijet events 

lA.eiii 

0.31 

Herwig-i--i- dijet events 

15.9/22 

0.82 

Data with 110 < mp < 140 GeV and 40 < mj2 < 60 GeV 

12.1/11 

0.79 

Data with 40 < nij < 60 GeV for both jets 

19.8/13 

0.56 

Data with 110 < nij < 140 GeV for both jets 

5.0/6 

0.91 


7 Systematic uncertainties 


The uncertainty on the background expectation is determined by the fitting procedure, which assumes a 
smoothly falling m,/ distribution. Possible uncertainties due to the background model were assessed by 
investigating several alternative families of parametrisations, and by considering signal plus background 
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Figure 3: Fits of the background model to the dijet mass (mjj) distributions in (a) Pythia 8 and (b) Herwig++ simu¬ 
lated background events that have passed all event selection and tagging requirements. The events are reweighted in 
both cases to correctly reproduce the leading-jet px distribution for untagged events, and the simulated data samples 
were scaled to correspond to a luminosity of 20.3 fb"^. The significance shown in the inset for each bin is calculated 
using the statistical errors of the simulated data. 




Figure 4: Fits of the background model to the dijet mass (nijj) distributions in data events (a) before boson tagging, 
and (b) where both jets pass all tagging requirements except for the nij requirement, and instead satisfy 40 < nij < 
60 GeV. 


fits of the chosen function to simulations of the dominant background as well as sidebands and control 
regions of data in which a signal contribution is expected to be negligible. These effects were estimated 
to be no more than 25% of the statistical uncertainty at any mass in the search region. The effect of 
the uncertainty on the trigger efficiency, the variations of the selection efficiencies as a function of the 
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kinematic properties of the background, and the composition of the background were also studied and 
were found to be well-covered by the uncertainties from the fit. 

Systematic uncertainties on the shape of the nijj distribution and the normalisation of the W and Grs 
signal are expressed as nuisance parameters with specified probability distribution functions (pdfs). The 
overall normalisation is a product of scale factors, each corresponding to an identified nuisance para¬ 
meter. If the shape is affected by a given nuisance parameter, the systematic change is included when the 
signal distribution is generated. If the nuisance parameter does not affect the shape, but only affects the 
normalisation, the distribution is simply scaled. 

The jet pj scale apj is defined as a multiplicative factor to the jet pj in simulation, pj = apj p^^- 
Following the technique used in ref. [47], the systematic uncertainty on apx is assessed by applying 
the jet reconstruction and filtering algorithms to inner-detector track constituents, which are treated as 
massless, and matching these track jets to the calorimeter jets. The ratio of the matched track jet’s pj to 
the calorimeter jet’s px as a function of several kinematic variables is compared in simulation and data 
and found to be consistent within 2%. Hence, the pdf used for apx is a Gaussian with a mean of one and 
a standard deviation of 0.02. Similar methods are used to determine the pdfs for the scale uncertainties in 
jet mass mj and momentum balance -^fy. 

Mismodelling the jet pj resolution can change the reconstructed width of a diboson resonance. The jet 
px resolution in the simulation is 5% in this kinematic region, and a 20% systematic uncertainty on this 
resolution is implemented by applying a multiplicative smearing factor, rp, to the px of each reconstructed 
jet with a mean value of unity and a width \o-r^\. The nuisance parameter represents the uncertainty 
in the px resolution of reconstructed jets and is assumed to have a Gaussian pdf with a mean of zero and 
a standard deviation 0.05 x 

The pdfs for the uncertainties in jet-mass resolution and momentum-balance resolution, listed in table 4, 
are similarly constructed. 

The ntrk variable is not modelled sufficiently well in simulation [50], so it is necessary to apply a scale 
factor to the simulated signal to correct the selection efficiency of the ntrk requirement. A scale factor of 
0.90 ± 0.08 is derived from the ratio of the selection efficiency of this cut in a data control region enriched 
with VT/Z -I- jets events, where a high-px IT or Z boson decays hadronically, to the selection efficiency 
of this cut in simulation. The data control region is defined by selecting events in the kinematic range 
where the jet trigger used in the search is fully efficient, and where only the leading-px jet passes the 
tagging requirement on Fits to the jet mass spectrum of the leading-px jet determine the number of 
hadronically decaying W and Z bosons reconstructed as a single jet that pass the ntrk requirement as a 
function of the selection criteria on ntrk- The dominant uncertainty on these yields is the mismodelling of 
the jet mass spectrum for non-IT or non-Z jets, and is evaluated by comparing the yields obtained when 
using two different background models in the fit. The resulting scale factor is 0.90 ± 0.08. Since the ntrk 
requirement is applied twice per event in the selections used in the search, a scale factor of 0.8 is applied 
per selected signal event, with an associated uncertainty of 20%. 

A 5% uncertainty on the signal efficiency due to uncertainties on the parton shower and hadronisation 
model is also included. The uncertainty is estimated by comparing the selection efficiencies obtained in 
simulated signal samples generated and showered with Pythia 8 to the selection efficiencies obtained in 
samples generated with Pythia 8 and showered with Herwig-i-i-. An additional 3.5% uncertainty on the 
signal acceptance due to uncertainties on the PDFs is considered. This uncertainty is estimated according 
to the PDF4LHC recommendations [51]. 
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Table 4 summarises the systematic uncertainties affecting the signal shape and the pdf constraining the 
associated nuisance parameter. The largest uncertainty on the shape of the reconstructed signal is due to 
the jet pt scale and resolution; the uncertainty in the scale introduces an uncertainty on the scale of the 
mass of the reconstructed resonant signal, and the resolution introduces an uncertainty of the width. The 
jet mass scale uncertainty also has an effect on the scale of the reconstructed mass, but this effect is less 
significant. Table 5 summarises the systematic uncertainties affecting the signal normalisation. The jet 
mass scale uncertainty affects both the shape and the normalisation. 

Table 4; Summary of the systematic uncertainties affecting the shape of the signal dijet mass distribution and their 
corresponding models. cr) in the table denotes a Gaussian distribution for the variable x with mean p and 

standard deviation cr. 


Source 

Uncerfainfy 

Consfraining pdf 

Jel pt scale 

2% 

G(apT| 1,0.02) 

Jel pt resolution 

20% 

G((r,^|0,0.05x Vl.22 - U) 

Jel mass scale 

3% 

G(an,|l,0.03) 


Table 5; Summary of the systematic uncertainties affecting the signal normalisation and their impact on the signal. 


Source 

Normalisafion uncertainly 

Efficiency of fhe Irack-mulliplicily cul 

20.0% 

Jel mass scale 

5.0% 

Jef mass resolution 

5.5% 

Subjel momenlum-balance scale 

3.5% 

Subjef momenlum-balance resolution 

2.0% 

Parlon shower model 

5.0% 

Parlon dislribufion funcfions 

3.5% 

Luminosily 

2.8% 


8 Results 

8.1 Background fit to data 

The fitting procedure is applied to the data after ITZ, WW and ZZ selection, and the results are shown in 
figure 5. In fhis figure, fhe fitfed background functions, labelled “background model,” are again integrated 
over fhe same bins used fo display fhe dafa, and fhe size of fhe shaded band refiecfs fhe uncerfainfies on 
fhe paramefers propagafed fo show fhe uncerfainfy on fhe expecfafion from fhe fif. Figure 5 also displays 
fhe filled dijel mass dislribufion of evenls passing any of fhe Ihree lagging seleclions. The lower panels in 
fhe figure show fhe significance of fhe difference belween dafa and fhe expecfafion in each bin. Table 6 
gives fhe tiffed values of fhe paramefers for fhe dafa selecled before fagging, displayed in figure 4(a), and 
afler fhe VTZ, WW and ZZ seleclions, as well as fhe number of evenfs observed. 
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Figure 5; Background-only fits to the dijet mass {mjj) distributions in data (a) after tagging with the WZ selection, 
(b) after tagging with the WW selection, (c) after tagging with the ZZ selection, and (d) for events passing any of 
the three tagging selections. The significance shown in the inset for each bin is the difference between the data 
and the fit in units of the uncertainty on this difference. The significance with respect to the maximum-likelihood 
expectation is displayed in red, and the significance when taking the uncertainties on the fit parameters into account 
is shown in blue. The spectra in the three signal regions are compared to the signals expected for an EGM W with 
mty/ = 1.5, 2.0, or 2.5 TeV or to an RS graviton with — 1.5 or 2.0 TeV. 


The dijet mass distributions after all three tagging selections are well-described by the background model 
over the entire mass range explored, with the exception of a few bins near nijj = 2 TeV which contain 
more events than predicted by the background model. Approximately 20% of the events selected by 
either the WW, WZ, or ZZ selection are shared among all three signal regions. The fraction of events 
common to the WZ and the WW or the WZ and the ZZ selections are 49% and 43% respectively. After 
requiring that mjj > 1.75 TeV, 5 out of 25 events are common to all three signal regions. The statistical 
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Table 6: Number of observed events and parameters from the background-only fits to the dijet mass spectrum for 
each tagging selection. The parameter ^ is a constant chosen after the fit to minimise the correlation between the 
fitted parameters p 2 and ■ 


Parameter 

Before tagging 

ITZ 

WIT 

ZZ 


4.3 

3.8 

4.2 

4.5 

Pi 

30.95 ± 0.03 

31.0 ± 1.4 

32.5 ± 1.5 

39.5 ± 2.0 

P3 

-5.54 ±0.03 

-9.1 ± 1.5 

-9.4 ± 1.6 

-9.5 ±2.3 

Observed events 

1335762 

604 

425 

333 


interpretation of these dijet mass distributions is discussed in the following section. 


8.2 Statistical analysis 

A frequentist analysis is used to interpret the data. For each of the two benchmark models under test, the 
parameter of interest in the statistical analysis is the signal strength, p, defined as a scale factor on the 
total number of signal events predicted by the model. Thus, the background-only hypothesis corresponds 
to ju = 0, and the hypothesis of a signal-plus-background model corresponds to p = 1. A test statistic, 
A{p), based on the profile likelihood rafio [52] is used fo tesf fhese models. The fesf stafisfic is designed to 
extract the information on p from a maximum-likelihood fit of the signal-plus-background model to the 
data. 

The likelihood model for the observation is 

X - n Ppoi,{nijni,p) X G(apT) x G(a^) x x N{d) (3) 

i 

where Ppoisin‘^^i,Jn'^xp) fhe Poisson probability to observe events if events are expected, G(ap.j.), 
GffTm), and GCcr^^), are the pdfs of the nuisance parameters modelling the systematic uncertainties related 
to the shape of the signal, and Af is a log-normal distribution for the nuisance parameters, 6, modelling 
the systematic uncertainty on the signal normalisation. The expected number of events is the bin-wise 
sum of the events expected for the signal and background: ngxp = ngig -i- nbg. The number of expected 
background events in dijet mass bin i, is obtained by integrating dnidx obtained from eqn. (1) over 
that bin. Thus nbg is a function of the dijet background parameters pi,p 2 ,P 3 - The number of expected 
signal events, ngig, is evaluated based on MC simulation assuming the cross section of the model under 
test multiplied by the signal strength and including the effects of the systematic uncertainties described in 
section 7. The expected number of signal events is a function of p and the nuisance parameters modelling 
the systematic uncertainties on the signal. 

The compatibility of the data with the background-only expectation is quantified in terms of the local po, 
defined as the probability of the background-only model to produce an excess at least as large as the one 
observed and quantified with an ensemble of 500,000 background-only pseudo-experiments, while the 
global probability of an excess with a given local po being the most significant excess to be observed any¬ 
where in the search region is quantified with 100,000 background-only pseudo-experiments that take into 
account the mass ranges and overlapping event samples for the three channels. The largest discrepancies. 
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in the region around 2 TeV in figures 5(a), 5(b) and 5(c), lead to small po values near that mass. The 
smallest local po values in the VTZ, WW, and ZZ channels correspond to significances of 3.4 cr, 2.6 cr, 
and 2.9 cr respectively. Considering the entire mass range of the search (1.3-3.0 TeV) in each of the three 
search channels, the global significance of the discrepancy in the WZ channel is 2.5 cr. 

Exclusion limits at the 95% confidence level are set following the CLs prescription [53]. 


8.3 Exclusion limits on new diboson resonances 

Limits on the production cross section times branching ratio of massive resonances are set in each diboson 
channel as a function of the resonance mass using the EGM W as a benchmark for the WZ selection, 
and the bulk Grs model for the WW and ZZ selections. In most of the mass range, the observed limit is 
somewhat better than the expected limit, but in the region near 2 TeV the excess of events in the data leads 
to observed limits which are weaker than expected. Eigure 6(a) shows the observed 95% CL upper limits 
on the cross section times branching ratio on the EGM W' WZ hypotheses as a function of the W' mass. 
EGM W' —> WZ for masses between 1.3 and 1.5 TeV are excluded at 95% CL. Eigures 6(b) and 6(c) show 
the observed 95% CL upper limits on the cross section times branching ratio for the bulk Grs ^ WW and 
ZZ, respectively. The cross section times branching ratio for excited graviton production with the model 
parameters described in section 3 is too low to be excluded with the sensitivity of this measurement. 


9 Conclusions 

A search has been performed for massive particles decaying to WW, WZ, or ZZ using 20.3 fb“^ of 
= 8 TeV pp collision data collected at the LHC by ATLAS in 2012. This is the first ATLAS search for 
resonant diboson production in a fully hadronic final state and strongly relies on the suppression of the 
dijet background with a substructure-based jet grooming and boson tagging procedure. The boson tagging 
selection includes different jet mass criteria to identify W and Z boson candidates and thus produces three 
overlapping sets of selected events for the searches in the WW, WZ, and ZZ decay channels. The most 
significant discrepancy with the background-only model occurs around 2 TeV in the WZ channel with a 
local significance of 3.4 cr and a global significance, taking the entire mass range of the search in all three 
channels into account, of 2.5 cr. 

Upper limits on the production cross section times branching ratio of massive resonances are set in each 
diboson channel as a function of the resonance mass, using an EGM W' WZ as a benchmark for 
the WZ channel, and an excited bulk graviton Grs to represent resonances decaying to WW and ZZ. A 
IT'with EGM couplings and mass between 1.3 and 1.5 TeV is excluded at 95% CL. 
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